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We investigate the dynamics of the generalized ACDM model, which the A term is running with 
the cosmological time. The A(t) term emerges from the covariant theory of the scalar field cj> with 
the self-interacting potential V((p). On the example of the model A(t) = Ab are + we show the 
existence of a mechanism of the modification of the scaling law for energy density of dark matter: 

Pdm oc a _3 +M*)_ yy e discuss the evolution of A(t) term and pointed out that during the cosmic 
evolution there is a long phase in which this term is approximately constant. This effect justifies 
Alcaniz and Lima’s approach to A (H) cosmologies. We also present the statistical analysis of 
both the A(t)CDM model with dark energy and decaying dark matter and the ACDM standard 
cosmological model. We divide the observational data into two groups: low z data (SNIa, BAO, 

H(z) and AP test) and high a data (Planck, WP and lensing). While for the former we find the best 
fit value of the parameter A is positive (A = 0.0338, energy transfer is from the dark energy to dark 
matter sector), for the latter we find that A is —0.0199 which is an evidence that the energy transfer 
is from decaying dark matter. This disagreement correlates with estimated values of Ho (67.77 
km/(s Mpc) and 65.62 km/(s Mpc) respectively). The decaying dark matter causes to lowering a 
mass of dark matter particles which are lighter than CDM particles and remain relativistic. The 
rate of the process of decaying matter is estimated. We show that in the models of decaying dark 
matter, the cosmological constant problem disappears naturally. The model with decaying dark 
matter possesses one parameter more but in light of the AIC it is better than the ACDM standard 
cosmological model. 
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I. INTRODUCTION 


In cosmology, the standard cosmological model (ACDM model) is an effective theory which well describes the 
current Universe in the accelerating phase of the expansion. All the astronomical observations of supernovae SNIa 
and measurements of CMB favor this model over the alternatives but we are still looking for theoretical models to 
dethrone the ACDM model. 


On the other hand the ACDM model has serious problems like the cosmological constant problem or the coincidence 
problem which are open and waiting for a solution. Among different propositions, it is an idea of introducing the 
running cosmological term |T]. The most popular way of introducing a dynamical form of the cosmological term is a 
parametrization by the scalar field, i.e. A = A(</>) or the Ricci scalar, i.e. A = A(R), where R is the Ricci scalar. In the 
both mentioned cases, the covariance of field equation is not violated and A = A(t) relation emerges from covariant 
theories i.e. A = A(</>(f)). For example in paper pQ, the relation A = Abare + jt is emerging from the theory of the 
scalar field with the self-interacting potential It is interesting that this type of a A(f) relation is supported by 

the non-critical string theory consideration [2 . 

Two elements appear in the ACDM model, namely dark matter and dark energy. The main aim of observational 
cosmology is to constrain the density parameters for dark energy as well as dark matter. In the testing of the ACDM 
model, the idea of dark energy is usually separated from the dark matter problem, The latter is considered as the 
explanation of flat galactic curves. Of course the conception of dark matter is also needed for the consistency of the 
model of cosmological structures but the hypothesis of dark energy and dark matter should be tested not as a isolated 
hypothesis. 

In this paper, we explore the A(f)CDM model with A(f) = Abare + f t , where t is the cosmological time for which 
we know an exact solution [l] . This enables us to show the nontrivial interactions between the sectors of dark matter 
and dark energy. It would be demonstrated that the model, which is under consideration, constitutes the special case 
of models with the interaction j] term Q = — ■ We will be demonstrated that the time dependence of the A term 

is responsible for the modification of the standard scaling law of dark matter pdm = Pdm.oo 3 ; where a is the scale 
factor [I] . Wang and Meng [3j developed a phenomenological approach which is based on the modified matter scaling 
relation p m = p m ,oa _3+<5 j where S is the parameter which measures a deviation from the standard case of cold dark 
matter (CDM). 

The both effect of the decaying A term and the modification of the scaling relation are strictly related in our model. 
One can obtain that CDM particles dilute more slowly in comparison to the standard relation p m oc a~ 3 in this model. 
In this context, the interesting approach was developed by Alcaniz and Lima 110. The coupling parameter S is also 
a subject of testing using astronomical data 


In this paper, due to it is known the exact solutions of our model it is possible to check how it works the model 
and one can strictly constrain the model parameters [T]. 

We estimate the value of A (t): pdm = Pdm,o« -3+A ^ where pdm is energy density of dark matter. In our methodology 
we consider the astronomical data as two separate data sets with respect to redshift 2 . The first group of low redshift 
consists of SNIa, BAO, H{z) and the AP test, and the second group is composed of Planck, WP and lensing data. 
The case of decaying dark matter effect is only relevant for high redshift data. This process takes place in the early 
Universe. 


We also analyze the model under considerations in details. In this analysis the model with A (t) = Ab are + % is our 
case study. For this model we show the terms A(f), S(t) are slow-changing with respect to the cosmological time and 
it justified to treat them as constants. 


The organization of the text is following. In Section III we present the model with A(f) = Abare + %- and its 
interpretation in scalar field cosmology with and witlrthe potential V(<f>). In Section III it is demonstrated how 
A(f)CDM cosmologies can be interpolated as interacting cosmologies with the interacting term Q = — d ^p and how 
they solve the coincidence problem (9|. In Section IV we present some results of the statistical estimations of the 
model parameters obtained from some astronomical data. Finally the conclusion are presented in Section [V] 


II. A(t)CDM COSMOLOGY WITH A — Abare "P , 2. 

Let us consider about the flat cosmological model with homogeneity and isotropy (the Robertson-Walker symmetry). 
The source of gravity is in the time dependent cosmological term and matter is in the form of a perfect fluid with 
energy density p m = p m (t ), where t is the cosmological time. The cosmic evolution is determined by the Einstein 
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equations which admit the Friedmann first integral in the form 

a 2 

3 H 2 (t) = Pm(t) + Abare + ) (1) 

where H{t) = dl ° [ f a is the Hubble function and a(t) is the scale factor and a 2 € M is a real dimensionless parameter. 
The sign of a 2 depends of the type of particle and the distribution of their energy. In the generic case the Breit-Wigner 
distribution gives rise the negative sign of a 2 unj. Note that this parametrization is distinguished by a dimensional 
analysis because a dimension of H 2 should coincide with a dimension of a time dependent part of A (t). 

It is assumed that the energy-momentum tensor for all fluids in the form of perfect fluid satisfies the conservation 
condition 


Tg = 0 , ( 2 ) 

where T a P = Tgf + A(t)g a P. The consequence of this relation is that 

p m + 3Hp m = - (3) 

at 

The cosmic evolution is governed by the second order acceleration equation 

H = -H 2 -^(p eS + 3p eB ) 1 (4) 

where p eB and p eB are effective energy density of all fluids and pressure respectively. In the model under the consid¬ 
eration we have 


Pe ff — Pm T PA i 


( 5 ) 


Pe ff — Pm Pa, 


( 6 ) 


where p m = 0, p\ = Ab are + pr and a 2 is a real number. 


For this case the exact solution of (JlJ) and (|3| for the Hubble parameter h = can be obtained in terms of Bessel 
functions 


T ( 3y/ClA,oHo + 

10 n , _ ln ~ l l 2 1 

h (t) = „„ + \/-fiA,o-7 — , - 

3 H 0 t v r f 3y/n A ,oH 0 i 


( 7 ) 


A ol 2 T 2 

where H 0 is the present value of the Hubble constant, Ua,o = -ypr, bI a ,o = 3H a , T 0 is the present age of the Universe 

T 0 = dt and n = 1 + i s index of the Bessel function. From |t|, the expression for the scale 

factor can be obtained in the simple form 


a(t) = C 2 

The inverse formula for t(a ) is given by 

e 

t(a) = ~ 


4.», 


~i —l 


y/2ni T ‘ 


3i^n K0 H 0 V 3 y/nIpH 0 C 2 


where S n (x ) is the Riccati-Bessel function S n (x) = y/^J n +k(x). 

Finally the exact formula for total mass p m {i) = pdm{t) + Pb(t) is given in the form 


/ 


Pm = ~3 Hq 


£2a,o 


Qo:,oTq 

i 2 


( 


1-2 n 
3 H 0 t 


L- 


y/n. 


A,0" 


/ 3^/n Ai0 H 0 \ \ 


| 3 ^a.o Ho t 


/ ) 


( 8 ) 


( 9 ) 


( 10 ) 
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Pdm(log(a)) 



FIG. 1. A diagram of the evolution of pdm(log(a)). The top thick line represents the evolution of pdm(log(a)) for Ho = 67.64 
km/(s Mpc), O m , 0 = 0.3696 and fl a 2 = 0.0392. The bottom thick line represents the evolution of Pdm(log(a)) for Ho = 68.13 
km/(s Mpc), f2 m ,o = 0.2646 and Q . a 2 = —0.0258. The medium line represents the best fit (see Table|I|. The gray region is the 
2cr uncertainties. We assumed 87 iG = 1 and we choose km J /(100 2 Mpc 2 s 2 ) as a unit of Pdm(log(a)). 

PdeOog(a)) 



FIG. 2. A diagram of the evolution of pdeQog(a)). The bottom thick line represents the evolution of pd e (log(a)) for Ho = 67.64 
km/(s Mpc), n m , 0 = 0.3696 and Q . a 2 = 0.0392. The top thick line represents the evolution of pd e (log(a)) for Ho = 68.13 km/(s 
Mpc), fi m ,o = 0.2646 and Q a 2 = —0.0258. The medium line represents the best fit (see Tabic |T|. The gray region is the 2 ct 
uncertainties. We assumed 8nG = 1 and we choose km 2 /(100 2 Mpc 2 s 2 ) as a unit of pd e (log(a)). 


,0 

of 1 


2 
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The diagram of pdm, Pde and P p^°^a)j as a f unc ti° n of log a obtained for low 2 data is presented in Fig 
Note that at the present epoch (log(a) = 0) both energy densities of dark matter and dark energy are of the same 
order (Fig. [2]) . 

While the relation A = A(i) violates the covariance of the general relativity Lagrangian, it can be simply demon¬ 
strated that such a relation can emerge from the covariant theory of scalar fields. This relation reveals as a time- 
dependent first integral of equation of motion. 

The action of general relativity for a perfect fluid has the following form 


s= + £m)d A x, 


(ii) 


where R is the Ricci scalar, C m = — p (\ + / ^^-dp^j [TT] and g is the metric tensor. The signature of g^ is chosen 
as -)■ 

For the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric without the curvature the Ricci scalar is expressed 
by R = + The Einstein equations are consequence of the variation of the Langragian C. The Einstein 

equations for dust and a minimal coupling scalar field are the following 

3 H 2 = p (12) 


and 


1 

6 


(p + 3p). 


a 


(13) 
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PJmHog(g)) 
Pde dog (a)) 



FIG. 3. A diagram of the evolution of ) ■ The top thick line represents the evolution of ^^1°^)) f° r = 67.64 

km/(s Mpc), fi m ,o = 0.3696 and fl a 2 = 0.0392. The bottom thick line represents the evolution of for Ho = 68.13 

km/(s Mpc), n m>0 = 0.2646 and fl a 2 = —0.0258. The medium line represents the best fit (see Tablejlj). The gray region is the 
2<j uncertainties. 


9 2 

We assume that p = p m + Pde and an interaction Q = between p m and pde and we assume the state equation for 

12 

p m as p m = 0 and for p^ e as p^e = —pde- Now we choice a new parametrization of pde as pd e = 2 + V” (</>), where 
cf> = Then (12) and (13) have the following form 


3 H 2 — p m + — + V (</>) 


and 


- = -\(p-p-2V 

a 6 V 

The conservation equation is in the form 

p + 3H (p + p) = p m + p de + 3 Hp m = 0. 

9 2 

Because there is an interaction Q = between p m and p^ e then the conservation equation gives two equations 

• 

Pm T d-Ti Pm — ,0 

t 6 


(14) 

(15) 

(16) 

(17) 


and 


Pde — 


dVU) ■ 2a 2 


(18) 


The first integral of (18) is ^- + V((j>) = jr + Abare- The potential V((j>) can be of any form. So if we want to find an 
exact solution for </>(fJwe should first assume the form of the potential V((j>). 


III. HOW A(i)CDM MODEL MODIFIES THE SCALING RELATION FOR DARK MATTER 

The existence of dark matter in the Universe is motivated by modern astrophysical observations as well as cosmo¬ 
logical observations. From observations of rotation curves of spiral galaxies, masses of infracluster gas, gravitational 
lensing of clusters of galaxies to cosmological observations of the cosmic microwave background anisotropy and large 
scale structures we obtain strong evidences of dark matter. 

Because models of nucleosynthesis of the early Universe are strongly restricted the fraction of baryons, we conclude 
that the nature of dark matter cannot be baryonic matter. On the other hand we imagine that particles of dark matter 
form a part of standard model (SM) of particles physics. There are many candidates for particles of dark matter 
e.g. WIMPs. Lately sterile neutrinos have been also postulated in this context mm- The interesting approach 
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is a search of photon emission from the decay or the annihilation of dark matter particles through the astrophysical 
observations of X-ray regions |141ll6] . For example the radiatively decaying dark matter particles as sterile neutrinos 
have been searched using X-ray observations US- 

Let us consider the ACDM model which describes a homogeneous and isotropic universe which consists of baryonic 
and dark matter and dark energy. Let us assume an interaction in the dark sectors. Then the conservation equations 
have the following form 


Pb + 3Hp h = 0, 
Pdm T 3HPdm — Q , 

Pde + 3 H pd e — 


(19) 

( 20 ) 
( 21 ) 


where pb is baryonic matter density, pd m is dark matter density and pd e is dark energy density {5j. Q describes the 
interaction in the dark sector. 

Let Pm = Pb + Pdm then (19) and (20) give 


P n 


3Hp m = Q. 


( 22 ) 


For model with A(f) = Ab are + tt the conservation equation has the form p m + 3 Hp m = — . 


So this model can 


be interpreted as the special case of model with the interacting in the dark sectors. In this model Q = — dl ^' > = 

n 2 rj _ 2 

Equation (221 for Q = -p- can be rewritten as p m + 3 H p m = Hp m t3Hp or 


dfm _ da ( 3 2a 2 \ 

p m a V t 3 Hp m ) ' 


(23) 


The solution of equation (23) is 


Pm — Pm,0& 


-3+<5(t) 


(24) 


where S = J S(t)d\oga, where 5(t) = t3H f^ p ^ ■ Q can be written as Q = S(t)Hp m . The evolution of <5(log(a)) 

and <5(log(a)) is presented in Fig. [ij and [ 5 ] One can observe that S(t) and 6(t) is constant since the initial singularity 
to the present epoch. 


<S(log(a)) 



log(a) 


FIG. 4. A diagram of the evolution of <5(log(a)). The top thick line represents the evolution of <5(log(a)) for Ho = 68.67 km/(s 
Mpc), O m , 0 = 0.3311 and Sl a 2 = 0.1034. The bottom thick line represents the evolution of ci(log(a)) for Ho = 67.34 km/(s 
Mpc), n m , 0 = 0.2975 and 2 = —0.0626. The medium line represents the best fit (see Table [T|. The gray region is the 20 - 
uncertainties. 


If S(t) is a slowly changing function than S(t) = S(t) = 6 and (241 has the following form 


Pm — Pm,oa 


— 3+5 


(25) 


In this case Q = SHp m . 

The early time approximation for S(t) is 


6(t) 


18 a 2 


4V1 + 3a 2 + 6a 2 + 4 


(26) 
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<5(log(a)) 



log(a) 


FIG. 5. A diagram of the evolution of 5(log(a)). The top thick line represents the evolution of 5(log(a)) for Ho = 68.67 km/(s 
Mpc), f2 m ,o = 0.3311 and = 0.1034. The bottom thick line represents the evolution of <5(log(a)) for Ho = 67.34 km/(s 
Mpc), S4 m ,o = 0.2975 and 12^ = —0.0626. The medium line represents the best fit (see Table [H). The gray region is the 20- 
uncertainties. Note that if pdm = 0 for a 2 < 0, i.e. whole dark matter decays then we have the ACDM model with baryonic 
matter. 


If 6(t) = 6 = const we can easily find that 

a = cio (27) 

Let us note that during the regime S = const both the time dependent parametrization and quadratic parametrization 
A (H) = Abare + I3H 2 coincide in Alcaniz and Lima approach [3]. 

We can rewrite pdm as 

Pdm = Pdm,oa -3+A ^ a \ (28) 

where A(t) = log n "n° a o _n^ lh ■ For the present epoch we can approximate A(f) as A (t) = A = const. So in the 
present epoch pdm = Pdm,oa -3+A . In the consequence, the Friedmann equation can be written as 3 H 2 = pbd^ 3 + 
PdmCi~ 3+x + Abare + ^2 ■ The evolution of A(log(a)) is presented in Kig. [gJ One can observe that A(f) is constant since 
the initial singularity to the present epoch. 


A(log(a)) 



log(a) 


FIG. 6. A diagram of the evolution of A(log(a)). The top thick line represents the evolution of A(log(a)) for Hq = 68.67 km/(s 
Mpc), 12 m ,o = 0.3311 and 12^ = 0.1034. The bottom thick line represents the evolution of A(log(a)) for Ho = 67.34 km/(s 
Mpc), 12 m ,o = 0.2975 and 12^ = —0.0626. The medium line represents the best fit (see Table [I]). The gray region is the 20- 
uncertainties. 


Following Amendola and others [181120] the mass of dark particles can be parametrized by the scale factor as 

m(a)=TOoexp J K(a , )d(loga'), (29) 


where mo is representing of mass of dark matter, k = . We consider the mass m(a) as an effective mass of 

particles in a comoving volume. 
















In Amendola et al. m the parameter k(o) is assumed as a constant. This simplifying assumption has physical justi¬ 
fication as it will be demonstrated i n th e further dynamical analysis of the model with decaying dark matter. Equation 
(29) can be simply obtained from (28) because e = a and m(a ) = a 3 p(a). Then A (a) = log * g j f a K(a')d(loga'). 
For illustration the rate of dark ma 


ter decaying process it would be useful to define the parameter fj 


/3 = 2^. (30) 

If A (a) = const then equation ( fip] ) has the equivalent form m(t) = moa^ exp~ log T , where /3 = 2^. 

Let consider the number of dark matter particles N(t) where t is the cosmological time. Than a half of the number 
of these particle N(t)/2 is reached at the moment of time /3f. 

For small value of 5 -C 3 the parameter /3 can be approximated by f3 = 2 . For example if we put A = 0.06 than 

the number of dark matter particles decreases to half after about 33.5 ■ 10 6 to starting at the time t 0 - 
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0 


-1 
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-2 



FIG. 7. A diagram of the switching phenomenon of a 2 . The bottom left line represents a 2 for the first fraction of particles and 
the top left line represents a 2 for the second fraction of particles before switching. The medium left line represents effective 
a 2 . After switching the bottom line represents a 2 of the first fraction which is equal zero. 

In Fig. [3] one can see that a quotient Pdm/ Pde decreases with the scale factor and remains of the same order for 
today (log a = 0). Therefore the coincidence problem is solved. 

If the dark matter sector consists of two fractions of particles than a 2 = a 2 s = af + a|, where a 2 is a value of 
a 2 for the first fraction of particles and a 2 is value of a 2 for the second fraction of particles. If particles of the first 
fraction is decayed than a 2 for this fraction should be equal zero and the value of a 2 ff is switched. Than a 2 ff = a 2 . 
This phenomenon is demonstrated in Fig. [7] 


IV. STATISTICAL ANALYSIS OF THE MODEL 


In this section, we present a statistical analysis of the model parameters such as H 0 , fldm.o and A. We are using 
the SNIa, BAO, CMB observations, measurements of H(z) for galaxies and the Alcock-Paczyhski test. 

We use the data from Union 2.1 which is the sample of 580 supernovae m- The likelihood function for SNIa is 

logisNia = -\[A- B 2 /C + \og{C/{2n))l (31) 

where A = ( p ohs - ^ th )C" 1 (At obs - ^ th ), B = C-\p ohs - p, th ), C = trC" 1 and C is a covariance matrix for SNIa. The 
distance modulus is p ohs = m — M (where m is the apparent magnitude and M is the absolute magnitude of SNIa) 
and /i th = 51og 10 D L + 25 (where the luminosity distance is D L = c(l + z) f* jjfzj)- 
We use Sloan Digital Sky Survey Release 7 (SDSS DR7) dataset at z = 0.275 j22], 6dF Galaxy Redshift Survey 
measurements at redshift z = 0.1 [23], the BOSS DR 9 measurements at z = 0.57 [23], and WiggleZ measurements at 
redshift z = 0.44, 0.60, 0.73 [25|. The likelihood function is given by 


log Lbao 



r s{Zd) \ r _i 

Dv(*)J 



r s (z d ) \ 


( 32 ) 


where r s (zd) is the sound horizon at the drag epoch [26) ■ 
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TABLE I. The best fit and errors for the estimated model for SNIa+BAO+I7(z)+AP test with Ho from the interval (65.0, 
71.0), f2dm,o from the interval (0.20,0.36) and Ao from the interval (—0.8,0.8). fib is assumed as 0.048468. The value of y 2 
for the best fit is equal 108.740 and the value of AIC is equal 114.74. In comparison with this model, the y 2 of the best fit of 
ACDM is equal 108.761 and AIC is equal 112.761. 


parameter best fit 

68 % CL 

95% CL 

Ho 

67.97 km/(s Mpc) 

+1.08 

-1.06 

+1.77 

-1.68 

^dm,0 

0.2650 

+0.0333 

+0.0561 

-0.0306 

-0.0488 

A 

0.0338 

+0.3832 

-0.3730 

+0.6406 

-0.6075 


The likelihood function for the Planck observations of cosmic microwave background (CMB) radiation [27], lensing 
from the Planck and low-1 1 polarization from the WMAP (WP) has the form 

log LcMB+lensing+WP = - \ (x‘ h - X° bs )C _1 (x th - X° bs ), (33) 

where C is the covariance matrix with the errors, x is a vector of the acoustic scale l a, the shift parameter R and 
f Ibh 2 where 


7iT [ Z dz' 

A ~VJ^) C J q H(z') 

R = y/n m , 0 Hg ^ ^ 7 ) 


here 2 * is the redshift of the photon-decoupling surface. 
The likelihood function for the Alcock-Paczynski test [5 


has the following form 


log Lap = 


(AP th ( Zi ) - AP obs (zi)Y 


(34) 

(35) 


(36) 


where AP(z) th = / 0 ‘ and AP(zj) obs are observational data [24 ] l30l - f37l . 

We are using some data of H{z) of different galaxies from 1381 - 140] and the likelihood function is 


log 



(37) 


We consider two cases. In the first case we estimate our model for low z astronomical observations such as SNIa, 
BAO, H(z) and the Alcock-Paczynski test. The final likelihood function for the first case is 

Ltot = -bsNia-bBAo Tap Lh (z) ■ (38) 

The second case is for high z and than we use Planck observations, lensing and WP. The value of Hq is chosen as 
Hq = 65.62 km/(s Mpc). 

We use our own code CosmoDarkBox in estimation of the model parameters. The code uses the Metropolis-Hastings 
algorithm mm and the dynamical system to obtain the likelihood function [261 43,. 

The results of statistical analysis for the first case are represented in Table 0 Figures [8] [9] and 1 10 | where it is shown 
the likelihood function with 68% and 95% confidence level projection on the (fidm,o> A) plane and the (U 0 , A) plane, 
respectively. 

For the second case the results of statistical analysis is shown in Table [IT] For the high z data we obtain that A is 
negative, which means that dark matter particles decay. 

We can use some information criteria in scientific practice to choose the best model. One of information criteria is 
the Akaike information criterion (AIC), which is given by 


AIC = —2 log L + 2d, 


(39) 
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FIG. 8. The likelihood function of two model parameters (Ha, Ao) with the marked 68% and 95% confidence levels for 
SNIa+BAO+iF(z)+AP test. The value of Hubble constant is estimated from the data as the best fit of value Hdm,o = 0.2650 
and then the diagram of likelihood function is obtained for this value. 


0.32 

0.30 

0.28 

I 

0.26 

0.24 

0.22 



FIG. 9. The likelihood function of two model parameters (fldm.o, Ao) with the marked 68% and 95% confidence levels for 
SNIa+BAO+Ff(s)+AP test. The value of the Hubble constant is estimated from the data as the best fit of value Ha = 67.97 
km/(s Mpc) and then the diagram of likelihood function is obtained for this value. 


where L is the maximum of the likelihood function and d is the number of model parameters. Model which is the 
best approximation to the truth from the set under consideration has the smallest value of the AIC quantity. It is 
convenient to evaluate the differences between the AIC quantities computed for the rest of models from our set and 
the AIC for the best one. Those differences AAIC 


AAIQ = AlCj - AIC m ,; n 


(40) 


are easy to interpret and allow a quick “strength of evidence” for the model considered with respect to the best one. 


V. CONCLUSIONS 

The main goal of the paper was to investigate in details the dynamic of the model with matter and the running 
cosmological constant term with respect to the cosmological time. It was assumed that baryonic matter satisfies the 
equation of the state for dust (i.e. is nonrelativistic). We were interested how the running A (t) influences on the 
scaling relation for energy density pdm- We have found the deviation from standard scaling a~ 3 for this relation. We 
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FIG. 10. The likelihood function of two model parameters (fidm.o, Ao) with the marked 68% and 95% confidence levels for 
Planck+lensing+WP. The value of the Hubble constant is equal 65.62 km/(s Mpc) and then the diagram of likelihood function 
is obtained for this value. 


TABLE II. The best fit and errors for the estimated model for Planck+lensing+WP for Ho = 65.62 s , with fldm.o from 
the interval (0.20,0.36) and Ao from the interval (—0.8, 0.8). fib is assumed as 0.048468. The value of y 2 for the best fit is 
equal 0.001 and the value of AIC is equal 4.001. The small value of y 2 is a consequence fitting of the two parameter model to 
the first peak of the spectrum of CMB. In comparison with this model, the value of y 2 of the best fit of ACDM for Ho = 0.6562 
is equal 13.929 and the value of AIC is 15.929. 


parameter best fit 68% CL 

95% CL 

Hdm,0 

0.2759 

+0.0068 

+0.0125 

-0.0072 

-0.0126 

A 

-0.0199 

+0.0057 

-0.0063 

+0.0102 

-0.0111 


explained the source of this deviation showing that pdm decreases more rapidly or slowly like a -3+<5 due to the energy 
transfer from dark matter to dark energy sector or in the opposite direction. The direction of the energy transfer 
crucially depends on the sign of a 2 constant in the model under consideration. 

The value of a 2 can be theoretically calculated in the quantum formalism developed by Urbanowski and Raczynska 
m- In their paper it was proposed a quantum mechanical effect which can be responsible for instability of charge 
particles which emit X or 7 rays. The sign of a 2 constant is obtained from the analysis of the decay law of unstable 
states of particles or its survival probability. In these calculations, the crucial role plays the Breit-Wigner distribution 
function which gives rise to a negative sign of the a 2 constant. For typical particles, decaying processes are describing 
through this distribution function. 

From the cosmological point of view it is interesting that fluctuations of energy of instantaneous states of particles 
of dark matter can be manifested as fluctuations of the velocity of particles [10] . In order this effect gives rise to 
the emission of the electromagnetic radiation from radio up to ultra-high frequencies. In the context of astrophysics 
important stays information can be obtained from the observation of X rays or 7 rays. From X ray CCD instruments, 
dark matter is searched in keV energy for looking for the non-baryonic X ray signature HI]. 

On the other hand the a 2 constant is a dimensionless model parameter which value can be estimated from some 
astronomical data. Our estimations favor the negative a 2 constant for high z, i.e. it is favored the decaying vacuum 
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of dark matter particles and the radiative nature of the energy transfer to dark energy sector. 

Sterile neutrinos give rise to a negative sign of the a 2 constant because they can be described in terms of the Breit- 
Wigner distribution function m ■ The negative sign of the a 2 constant offers a new insight into the cosmological 
constant problem because the running A is the growing function of the cosmological time with asymptotic Ab are at 
t —> oo. Therefore this problem disappears in a natural way. 

In our model the parameter a 2 cannot change a sign. On the other hand estimations showed the positive a 2 
for low redshift data and the negative one for high redshift data. One way of explanation of this discrepancy is to 
treat the parameter a 2 as some effective parameter describing two-component dark matter. For one component the 
contribution to a 2 is negative and for second one is positive. During the cosmic evolution the former decays entirely 
and the latter survives and there is the switch of the effective a 2 sign. 

In our paper we have also found the physical background of the approach developed by Alcaniz and Lima [1] in 
which the relation p^ m oc a _3+A , where A = const, plays an important role. Our observational analysis of the evolution 
this parameter during the cosmic evolution indicates that such an ansatz has a strongly physical justification. We 
showed using astronomical data that the model with decaying dark matter is favored on a 2cr level over the ACDM 
standard cosmological model. 

The good news are also coming from our consideration. In interacting cosmology the interacting term which is 
postulated in different physical forms is interpreted as a kind of nongravitional interactions in the dark sector. We 
suggest that this interaction has the radiation nature and can be rather interpreted following the Urbanowski and 
Raczynska idea as a possible emission of cosmic X and 7 rays by unstable particles [TO]. 

It is still an open discussion about the nature of dark matter: cold or warm dark matter [33]. Our results showed 
that in the model of dark matter decay dark matter particle being lighter than CDM particles. Therefore particles of 
warm dark matter remain relativistic longer during the cosmic evolution at the early universe. 

In the investigating dynamics of the interacting cosmology the corresponding dynamical systems, which are deter¬ 
mining the evolutional paths, are not closed until one specify the form of the interacting term Q. Usually this form 
is postulated as a specific function of the Hubble parameter, energy density of matter or scalar field or their time 
derivatives [451148] . The case study of our model with decaying dark matter favor the choice of the interacting term 
in the form Q oc Hp m . Zilioti et al. has recently demonstrated that the coincidence problem can be naturally solved 
in a large class of running A (H) vacuum cosmologies [49] . The case study of our model fully confirms this idea. 


Acknowledgements 

The work was supported by the grant NCN DEC-2013/09/B/ST2/03455. The authors thank prof. A. Borowiec 
and A. Krawiec for remarks and comments. 


[1] M. Szydlowski and A. Stachowski, (2015), arXiv: 1507.02114 [astro-ph.CO] 

[2] J. L. Lopez and D. V. Nanopoulos, Mod. Phys. Lett. All, 1 (1996), arXiv:hep-ph/9501293 [hep-ph] 

[3] P. Wang and X.-H. Meng, Class. Quant. Grav. 22, 283 (2005), arXiv:astro-ph/0408495 [astro-ph] 

[4] J. S. Alcaniz and J. A. S. Lima, Phys. Rev. D72, 063516 (2005), arXiv:astro-ph/0507372 [astro-ph] 

[5] L. L. Graef, F. E. M. Costa, and J. A. S. Lima, Phys. Lett. B728, 400 (2014), arXiv:1303.2075 [astro-ph.CO] 

[6] D. Bessada and O. D. Miranda, Phys. Rev. D88, 083530 (2013), arXiv:1310.8571 [gr-qc] 

[7] M. Szydlowski, Phys. Rev. D91, 123538 (2015), arXiv: 1502.04737 [astro-ph.CO] 

[8] R. S. Goncalves, G. C. Carvalho, and J. S. Alcaniz, (2015), arXiv:1507.01921 [astro-ph.CO] 

[9] S. del Campo, R. Herrera, and D. Pavon, Phys. Rev. D91, 123539 (2015), arXiv:1507.00187 [gr-qc] 

[10] K. Urbanowski and K. Raczynska, Phys. Lett. B731, 236 (2014), arXiv: 1303.6975 [astro-ph.HE] 

[11] O. Minazzoli and T. Harko, Phys. Rev. D86, 087502 (2012), arXiv: 1209.2754 [gr-qc] 

[12] A. Boyarsky, O. Ruchayskiy, and M. Shaposhnikov, Ann. Rev. Nucl. Part. Sci. 59, 191 (2009), arXiv:0901.0011 [hep-ph] 

[13] H. Motohashi, A. A. Starobinsky, and J. Yokoyama, Phys. Rev. Lett. 110, 121302 (2013), arXiv: 1203.6828 [astro-ph.CO] 

[14] K. Abazajian, G. M. Fuller, and W. H. Tucker, Astrophys. ,J. 562, 593 (2001), arXiv:astro-ph/0106002 [astro-ph] 

[15] A. Boyarsky, O. Ruchayskiy, D. Iakubovskyi, and J. Franse, Phys. Rev. Lett. 113, 251301 (2014), arXiv: 1402.4119 [astro- 
ph.CO] 

[16] A. D. Dolgov and S. H. Hansen, Astropart. Phys. 16, 339 (2002), arXiv:hep-ph/0009083 [hep-ph] 

[17] N. Sekiya, N. Y. Yamasaki, and K. Mitsuda, (2015), arXiv:1504.02826 [astro-ph.HE] 

[18] L. Amendola, G. Camargo Campos, and R. Rosenfeld, Phys. Rev. D75, 083506 (2007), arXiv:astro-ph/0610806 [astro-ph] 

[19] E. Majerotto, D. Sapone, and L. Amendola, (2004), arXiv:astro-ph/0410543 [astro-ph] 

[20] R. Rosenfeld, Phys. Lett. B624, 158 (2005), arXiv:astro-ph/0504121 [astro-ph] 

[21] N. Suzuki et al., Astrophys. J. 746, 85 (2012), arXiv:1105.3470 [astro-ph.CO] 




13 


[22] W. J. Percival et al. (SDSS), Mon. Not. Roy. Astron. Soc. 401, 2148 (2010), arXiv:0907.1660 [astro-ph.CO] 

[23] F. Beutler, C. Blake, M. Colless, D. H. Jones, L. Staveley-Smith, L. Campbell, Q. Parker, W. Saunders, and F. Watson, 
Mon. Not. Roy. Astron. Soc. 416, 3017 (2011), arXiv:1106.3366 [astro-ph.CO] 

[24] L. Anderson et al, Mon. Not. Roy. Astron. Soc. 427, 3435 (2013), arXiv: 1203.6594 [astro-ph.CO] 

[25] C. Blake et al., Mon. Not. Roy. Astron. Soc. 425, 405 (2012), arXiv: 1204.3674 [astro-ph.CO] 

[26] D. J. Eisenstein and W. Hu, Astrophys. J. 496, 605 (1998), arXiv:astro-ph/9709112 [astro-ph] 

[27] P. A. R. Ade et al. (Planck), Astron. Astrophys. 571, A16 (2014), arXiv: 1303.5076 [astro-ph.CO] 

[28] C. Alcock and B. Paczynski, Nature 281, 358 (1979) 

[29] M. Lopez-Corredoira, Astrophys. J. 781, 96 (2014), arXiv:1312.0003 [astro-ph.CO] 

[30] P. M. Sutter, G. Lavaux, B. D. Wandelt, and D. H. Weinberg, Astrophys. J. 761, 187 (2012), arXiv: 1208.1058 [astro- 
ph.CO] 

[31] C. Blake et al., Mon. Not. Roy. Astron. Soc. 418, 1725 (2011), arXiv: 1108.2637 [astro-ph.CO] 

[32] N. P. Ross et al., Mon. Not. Roy. Astron. Soc. 381, 573 (2007), arXiv:astro-ph/0612400 [astro-ph] 

[33] C. Marinoni and A. Buzzi, Nature 468, 539 (2010) 

[34] J. da Angela, P. J. Outram, and T. Shanks, Mon. Not. Roy. Astron. Soc. 361, 879 (2005), arXiv:astro-ph/0505469 
[astro-ph] 

[35] P. J. Outram, T. Shanks, B. J. Boyle, S. M. Croom, F. Hoyle, N. S. Loaring, L. Miller, and R.. J. Smith, Mon. Not. Roy. 
Astron. Soc. 348, 745 (2004), arXiv:astro-ph/0310873 [astro-ph] 

[36] I. Paris et al., Astron. Astrophys. 548, A66 (2012), arXiv:1210.5166 [astro-ph.CO] 

[37] D. P. Schneider et al. (SDSS), Astron. J. 139, 2360 (2010), arXiv: 1004.1167 [astro-ph.CO] 

[38] J. Simon, L. Verde, and R. Jimenez, Phys. Rev. D71, 123001 (2005), arXiv:astro-ph/0412269 [astro-ph] 

[39] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski, and S. A. Stanford, JCAP 1002, 008 (2010), arXiv:0907.3149 [astro- 
ph.CO] 

[40] M. Moresco et al., JCAP 1208, 006 (2012), arXiv: 1201.3609 [astro-ph.CO] 

[41] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, J. Chem. Phys. 21, 1087 (1953) 

[42] W. K. Hastings, Biometrika 57, 97 (1970) 

[43] W. Hu and N. Sugiyama, Astrophys. J. 471, 542 (1996), arXiv:astro-ph/9510117 [astro-ph] 

[44] A. Lapi and L. Danese, (2015), arXiv: 1508.02147 [astro-ph.CO] 

[45] L. Amendola, Phys. Rev. D62, 043511 (2000), arXiv:astro-ph/9908023 [astro-ph] 


[46] W. Zimdalil and D. Pavon, Phys. Lett. B521, 133 (2001), arXiv:astro-ph/0105479 


[47] L. P. Chimento, Phys. Rev. D81, 043525 (2010), arXiv:0911.5687 [astro-ph.CO] 



[48] M. Szydlowski, T. Stachowiak, and R. Wojtak, Phys. Rev. D73, 063516 (2006) 

arXiv:astro-ph/0511650 [astro-ph] 


[49] G. J. M. Zilioti, R. C. Santos, and J. A. S. Lima, (2015), arXiv: 1508.06344 [gr-qc] 



